Models of mesoscopic time-resolved current detection 
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Quantum transport in mesoscopic conductors is essentially governed by the laws of quantum me- 
chanics. One of the major open questions of quantum mechanics is what happens if non-commuting 
observables are measured simultaneously. Since current operators at different times do not com- 
mute, the high-frequency correlation functions of the current are realization of this fundamental 
quantum question. We formulate this problem in the context of measurements of finite-frequency 
current cumulants in a general quantum point contact, which are the subject to ongoing experimen- 
tal effort. To this end, we present two models of detectors that correspond to a weak time-resolved 
measurement of the electronic current in a mesoscopic junction. In both cases, the backaction of the 
detector leads to observable corrections to the current correlations functions involving the so-called 
noise susceptibilities. As a result, we propose a reinterpretation of environmental corrections to the 
finite-frequency cumulants as inevitable effect resulting from basic quantum mechanical principles. 
Finally we make concrete predictions for the temperature-, voltage-, and frequency-dependence of 
the third cumulant, which could be verified directly using current experimental techniques. 
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I. INTRODUCTION 

Since the advent of mesoscopic physics, the quantum 
mechanical properties of charge transport in conductors 
are under intensive research. Most manifestations to date 
are concerned with properties following from the quan- 
tum mechanical wave nature>i More recently the quan- 
tum mechanics of the electron spin has been investigated 
in electron transportfS. On the other hand, the aspect of 
measurement has attracted less attention so far in quan- 
tum transport. This is surprising in view of it's emi- 
nent importance for the fundamental difference to classi- 
cal physics.^ 

The electronic transport through mesoscopic junctions 
is quite well described by independent electrons.** Due to 
the Landauer formula, the conductance can be expressed 
in terms of the single-particle scattering matrix and the 
level occupations of the leads. The same scattering ma- 
trix can describe low-frequency current noise." Experi- 
mentally, the scattering description has been confirmed 
in the case of low frequency noised a while ago. 

It is remarkable, that the same scattering matrix also 
describes higher order correlation functions of current ~ 
leading to the so-called full counting statistics (FCS)i^ 
These results have been confirmed for the third cumu- 
lant of tunnel junctions only, so fari^i^ However, a proper 
account of the first experiments showed that the mea- 
sured third cumulant was governed by an environmental 
contribution^^ 

In the last years it has also been realized, that the 
measurement of the quantum noise depends in an essen- 
tial way on the detection schemed To understand this, 
one considers different Fourier transform of the current- 
current correlation function C{t,t') ~ {I{t)I{t')). It 
turns out that the unsymmetrized correlators C{uj) — 
J dTC{0, t) exp(±iwT) are related to absorption or emis- 
sion of an energy quantum huj, respectively, and can 
therefore be measured in a corresponding detector. On 



the other hand, a classical detector measures always the 
real combination C{io) + C{~uj), corresponding to the 
anti-commutator of the current operators. The high fre- 
quency noise has been also tested experimentallyjii^rJi 
The results are in full agreement with the scattering the- 
ory of non-interacting quasiparticles in mesoscopic con- 
ductors. 

From the above it is clear that a time-resolved detec- 
tion of the quantum mechanical current in a quantum 
transport process has to properly account for the quan- 
tum rules for measurements. This is not so much an 
experimental question - since experimentalists automat- 
ically obey the rules of quantum mechanics - rather, the 
proper treatment of a quantum measurement represents 
a challenge to theorists, who have to worry, which ex- 
pectation values are actually measurable. Hence, the ex- 
tension of full counting statistics to time-resolved (high- 
frequency) correlations requires first a proper definition 
in terms of projective measurement*? Qualitatively, one 
can distinguish two extremes. A continuous projective 
measurement would lead up in suppression of the dy- 
namics of the system - a phenomenon known as quan- 
tum Zeno effect jiS, To avoid the Zeno effect it is neces- 
sary to include the detector's degrees of freedom into the 
complete evolution. Some effects of the detector backac- 
tion have been addressed already in the literatureJ^"— 
However, instead of including a specific detector's dy- 
namics one can restrict the treatment to the system 
only by replacing the projection by so-called Kraus 
operators*22, These define a positive operator- valued mea- 
sure (POVM)^ and by virtue of Naimark's theorem the 
two methods are equivalents^ as long as initially the de- 
tector and the system are uncorrelated--'^ 

Once we have established the detection model we can 
again express current correlations of arbitrary order in 
terms of the single-particle scattering matrix. However, 
it turns out that a simple evaluation of these averages 
quickly leads to very complicated and cumbersome ex- 
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pressions. Nevertheless, the demanding task necessary to 
analyze the generating functional leads to interesting an- 
alytic results for the frequency-dependent cumulants and 
the full counting statistics.™"™ Below, we will provide a 
simplified framework, which allows to evaluate correla- 
tors of arbitrary order in an efficient way. 

An important recent development is that the third cu- 
mulant has been measured at high frequency.^'' These 
quite remarkable results show interesting quantum fea- 
tures at eV = fvjj. However, most of this feature 
can be explained by environmental effects due to a se- 
ries impedance. These result in a mixing of the third- 
order correlator with the noise susceptibilities due to 
the frequency-dependent backaction of an environmental 
impedance. After fitting the spurious contribution, the 
resulting third cumulant of the considered tunnel junc- 
tion shows no frequency dispersion at all. This is in per- 
fect accordance with the theory, which predicts in fact 
a dispersionless third cumulant in the case of a tunnel 
junctioui^ However, it is most likely just a question of 
time until the third cumulant for a mesoscopic conductor 
with non-opaque channel is measured. Hence, for these 
experiments it is vital to distinguish the environmental 
contribution from the quantum mechanical backaction of 
the detector. 

Our findings can be summarized in two main achieve- 
ments. On one hand, most of the article is devoted to 
the development of concrete models of mesoscopic detec- 
tors of high-frequency current-detection. The interplay 
between the quantum mechanical projection and backac- 
tion plays here a major role and we clariiy below, how 
this can or cannot be distinguished from the environ- 
mental effect. On the other hand, we present an efficient 
method to calculate higher order correlation functions of 
the current in mesoscopic junctions. This allows us to 
present several results for the third cumulant and the 
noise susceptibility in a unified and transparent manner. 
The method can be the basis for an efficient evaluation of 
fourth and higher order correlators, which is important 
for example in photon counting statistics 

In the main body of the article we are going to con- 
struct a model of a projective time-resolved measurement 
of the current in a mesoscopic conductor. We present 
two possibilities. One is to define a POVM, which corre- 
sponds to a weak measurement of the current in a quan- 
tum point contact. This phenomenological and very sim- 
ple approach has the appeal that the competition be- 
tween projection and weak measurement is qualitatively 
reflected. However, a concrete physical realization is dif- 
ficult to find. Another possibility is based on a phys- 
ical model of a detector, which we call quantum tape. 
Our detector will be a massless bosonic field described 
by charge density and current operators, parametrized 
by a one-dimensional variable. It is equivalent to the 
noninteracting Schwinger-Tomonaga-Luttinger model of 
a one-dimensional fermionic field. '^'^ Equivalently the de- 
tector can be thought of as a transmission line, which is 
described by the same Hamiltonian. This detector model 



is actually a massless version of a Josephson transmission 
line detector The field can be decomposed into left and 
right going components with a constant velocity. In this 
way, the spatial coordinate of the detector, which is cou- 
pled to the system only at one point in space, corresponds 
effectively to a time-coordinate. The information regis- 
tered by the detector is moved away from the system due 
to the internal dynamics of the detector. Finally a spa- 
tially resolved strong projection is applied to the detector 
at one instant of time, and the result can be translated 
into the temporal current profile. The projection can be 
made on one or both ends of the tape due to its bidirec- 
tional information transfer. The results always depend 
on the detection scheme. In the limit of weak coupling 
the detector's contribution is only a large Gaussian off- 
set noise, independent of the measured system. '^^ Accord- 
ingly, the influence of the backaction of the detector onto 
the system is also weak but always present. 

The secondary goal of the present article is to develop 
a compact method to calculate arbitrary cumulants of 
current operators in a junction with energy-independent 
transmission. To this end, we set up a current algebra, 
which is especially suited for non-equilibrium scattering 
problems. This becomes necessary, when considering dif- 
ferent orderings of higher-order correlators at finite fre- 
quency. A straightforward evaluation using the scatter- 
ing states is possible, but becomes in general very cum- 
bersome. Introducing an algebra for certain current op- 
erators of a quantum point contact greatly simplifies this 
calculation and will therefore be useful also in other con- 
text, hopefully. In the respective limits, our results are, 
of course, in agreement with existing results. 

The paper is organized as follows. We begin with the 
description of the mesoscopic junction in Sec. II, together 
with an overview of necessary mathematical tools. In Sec. 
Ill we present the POVM detection model and extend it 
to many detectors. In this model, we are able to derive 
also some higher cumulants of current in the limit of weak 
coupling between the detector and the system. In certain 
limits (weak transmission or refiection) one can find the 
complete PCS generating functional. Finally in Sec. IV, 
we present a strictly projective measurement model - the 
quantum tape, which gives results similar to the POVM. 
Many useful but lengthy mathematical details are moved 
to appendices. 



II. MESOSCOPIC JUNCTION 



The junction is defined as constriction in a two dimen- 
sional electron gas, which is narrow in the y-direction 
and relatively long in the x-direction. We assume no 
electron-electron interactions and no magnetic field. The 
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Hamiltonian of the system reads 



= ^ / dxdy -ipl{x,y)h'^Atpcr{x,y)/2Tn 



+eVe {x, y, 1)11)1 (x, y)'ip„ {x, y) 



(1) 



The fermionic operators satisfy standard anticommuta- 
tion relations 

{■ipa{x,y),ipa'{x',y')} = 0, 

{Mx, y),^l ix', y')} = 6,,,6ix - x')Siy ~ y'). (2) 



The external potential splits into three parts 
Ve{x,y,t)^V^{x) + Vy{y) + V{x,t), 



(3) 



where Vx is the scattering potential, Vy is the effect of the 
constriction and V is the time-dependent bias potential. 
The scattering potential is assumed to be nonzero only 
in a small interval around the center of the constriction 
while Vy grows to infinity with \y\ — oo. We have the 
total current operator in x direction 

m fdy (4) 

2m J 

^[4't{x,y)dx-4'a{x,y) - {dxiljU.x,y))'4'a{x,y)]- 

a 

In absence of a bias voltage, the Hamiltonian can be 
diagonalized in the space of scattering states described 
by fermionic operators i)k,n{x) with the channel quantum 
number including the electron spin n — (n, a). Far from 
the region 14 ^ 0, we have 



i^a{x,y) 



dk 



(5) 



Here we distinguish between the normalized transversal 
mode wavefunction ipn and the mode operator f/'fe.n, k is 
the longitudinal wave vector, f„ and in are reflection and 
transmission scattering amplitudes, respectively. Due to 
unitarity, transmission r„ = |i„p and reflection coeffi- 
cient Rn — l^^riP satisfy i?„ + — I. We assume that 
the junction is long enough to treat A: as a continuous 
parameter. The transverse mode index n is kept as dis- 
crete. The energy has the structure E — En + h^k"^ /2m, 
where En is for the transversal part. 



A. Expansion around Fermi level 

We are interested only in phenomena in a narrow part 
of the electron band around Fermi level Ep .— In particu- 
lar, we assume fc^T, eV, h\uj\ <^ Ep- We also neglect the 
exact structure of wavefunctions in the region Vx ^ Q and 
take only asymptotic states, like in ([5]). We denote Fermi 



wave numbers fc„ = yj2m{Ep — En)/h and Fermi veloc- 
ities Vn = hkn/m. Only modes with En < Ep contribute 
at zero temperature. 

We construct an extended Hilbert space consisting of 
left and right going states, L and R, respectively. The 
relation between operators in the reduced and standard 
space is 



'4^Lnix) 
■^Rn {x) 



fc<0 



dkdx'dv 
dkdx'dy 



■''^''-^H*n{y)M^',y), 



^ik(x' -x) I * 



ik>o 27r 

In fact, only k —kn and k kn play a role for L 
and R, respectively. The actual dynamics of states deep 
below or above the Fermi sea can be ignored and its only 
reminder will be some ultraviolet regularization or cutoff. 
The new operators satisfy anticommutation relations 

{'lpAri{x),i^Brh{x')} = 0, 

{V^Aftla;), -03^(2;')} = ^AB^ninK^ - x') (7) 

for A,B^ L,R. 

With help of the above stated approximations, we can 
write H in the new space as 



Ha + ^ QnVn 
n 

/ dxeV{x,t) 



i>L mm (0) + (0)^i« (0) + (8) 



^{n{^)'^Ln{x) + i'Rni^)i}Rni.^) 



where Hq is equal to 



dxihVn ■4^\n{x)dx'4>Ln{x) ~ i)'^p.n{x)dxi^Rn{x) 



(9) 

It is necessary to regularize the second term in Eq. ([8]) 
to define the transmission and the reflection coefficients, 
Tn — cos^(g„/fi.) and i?„ = sin^(q„/?i), respectively (see 
the discussion in the appendix A). 

The current operator ([5]) in the new basis is replaced 

by 

-^(^) XI ew«(V'k(2;)V'i?.ft(a;) - i^{f^{x)i)Ln{x)). (10) 

n 

As we will later see, it is disadvantegous to operate di- 
rectly with the field operators. To circumvent this it is 
possible to introduce bosonic operators (viz. quadratic 
forms of Fermion operators) and to develop a closed al- 
gebra for those. Since current operators are generally 
non-commuting, we have to cope with their algebra. To 
this end, we introduce the following auxiliary operators 

hn{Xn) = ^{i^\n{x)'ipLfXx) + i^m{-x)'pRn{-x)), 



^l7i (-^n ) 



-{■<pLni^)^Ln{x) - (-x) (-2;) ) , 



eVr, 



-^Lni^)^Rn{-x) +h.C. , 
■'^Lni^)'^Rn{~x) + h.C. . 



(11) 
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with Xn = x/vn- Note that Xn has the unit of time. It 
wih be more convenient for us than length units. The 
Hamiltonian ([8]) can be now written as 



H = Ho 



2g„/3fi(0)/e 



+ J2 [ dsVi±SVn,t)[ionis)±Iln{s)]. (12) 
n,± 

The current takes the form 

n 

(13) 

Using arguments similar to bosonization^S (see appendix 
A), we get the following commutation rule 

[ijn{s)Jkm{s')] = ie^5n,n5jkds5{s - s')/47r 

+ ^^if-jkie5nfhiin{s)5{s - s') (14) 

with ejki equal +1 for jkl — 123,231,312, —1 for jkl = 
321, 213, 132 and zero otherwise. The last useful commu- 
tator is 



[HoJ]n[s)] = -ihdsijn{s) 

and one can in principle write 



(15) 



(16) 



The current /q corresponds to the total charge changes 
in the leads. It is preserved due to absence of capacitive 
effects and behaves like a free bosonic field. Without the 
scattering term I3 in ([T^ only the first term in (|14p needs 
to be considered in the dynamics, since only Iq and Ii are 
dynamical variables. The system could be fully bosonized 
and all correlation functions can be calculated. A non- 
vanishing scattering retains some fermionic features due 
to the commutator [/i,/3]. As a consequence the fluctu- 
ations become non-Gaussian. 



B. Equilibrium averages 

Far from the junction and without bias potential, we 
can write the initial density matrix as 

p ^ cM-Ho/kBT)/TTeM-Ho/kBT). (17) 

The commutator gives the useful relation 

Trp/,^(s)i = Ti-pAijn{s + ih/kBT). (18) 

A similar relation has been introduced in Ref. IH. For 
convenience we rewrite the previous equation and the 
relations (|14l) in the frequency domain for the operators 



(19) 



From now on, a Greek argument will always denote 
the Fourier-transformed operators. The current algebra 
reads now 



[ijn{a)Jkrn{l3)] = e'^ SnrnSjkCtSia + f3)/2 

+ ^ if- jkl eSfirhlin (a + /3) • 



(20) 



Eq. (jl8p reproduces the fiuctuation-dissipation 
theorem^^ 



TTpIjniLo)A = e^BTTrpAIjf.iio). 



(21) 



With the above derived tools it is now straightforward to 
calculate equilibrium averages. As examples we obtain 



Tr pi jn{uj) = 0, 



(22) 



TTpI,^{a)Ik^{f3) = '-^^^S{a + f3){w{a) + a), 

Trp/,«(a)4^(/3)/,,(7) = ^fllM^l^Sia + /3 + 7) 
[u{l3){w{a) — w{j)) + w{a) — w{P) + w["f) -I- a — 7] , 

where we introduced 

u{uj) = cot'h{tiuj/2kBT) , (23) 
w{uj) — IjJu{lo) , 

Tr[kBT/hY 



w{s) = — lim Re— 



e^o sinh^ {nskBT/h + ie) 

We note that w(s — 0) is singular and all practical calcu- 
lations need a proper regularization around s = 0. This 
can for example be an ultraviolet cutoff of the order of 
the Fermi energy in the frequency domain. 



C. Heisenberg equations 

To generalize the previous method to a nonequilib- 
rium situation, it is necessary to solve the Heisenberg 
equations for the various current operators we have in- 
troduced. It is a major advantage of the current algebra, 
we have introduced, that this can be done exactly. From 
now on, we switch to Heisenberg picture, so all operators 
will be transformed as 



indtA{t) = [A{t),H(t)]+ihdtA{t). 



(24) 



Here the dt denotes explicit time dependence of the op- 
erator. Remarkable, the Heisenberg equations can be 
solved completely for the operators Ijn. The details are 
described in the appendix B and in the following we 
merely present the results. To distinguish between time 
and spatial coordinate (which has the same unit, when 
divided by Fermi velocity) we will denote by A{s,t) an 
operator at position s and time t. We assume an arbitrar- 
ily time-dependent bias potential V{x,t) = 9{—x)V{t) 
and neglect some regularization issues around x — and 
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X — >■ — oo, which can of course be properly handled us- 
ing the procedures mentioned previously. We find for the 
Heisenberg operators 



Ion{s,t) = Iofi{t + s), 

hn{s,t) = iifi{t + s) for s > 0, 

Iln{s,t) = (T„ - Rn)iinit + s) 



(25) 



'TTT^^^^i* + - WRnTnJ^it + s) for s < 



with 



= hni-s) cos$(s) - hn{s) sin$(s). (26) 



and 



$(s) 



dt eV{t)/h. 



(27) 



Similar expressions can be obtained for l2n{s,t) and 
Isfiisjt) and are presented in the appendix B. The big 
advantage of the above equations is that the Heisen- 
berg current operators are represented by linear combi- 
nations of equilibrium operators. Hence, their averages 
are straightforwardly obtained using the results of the 
previous subsection. 

We conclude this section by emphasizing the main re- 
sults we have obtained so far. The calculation of higher 
order cumulants of the current in a mesoscopic junc- 
tion using the standard scattering approach is quite a 
demanding task. This is true in particular at finite fre- 
quency, when operator ordering issues become even more 
prominent. The method, we have developed in this sec- 
tion builds in an elegant way on a very simple algebra of 
current operators (interestingly resembling the usual an- 
gular momentum algebra) . Our method allows to obtain 
in a relatively straightforward manner arbitrarily ordered 
current cumulants of (almost) any order. 



III. THE MEASUREMENT - POVM 

We now turn to the main topic of the article, how 
the current correlators at high frequencies can be mea- 
sured and how the measurement protocol itself influences 
the measurement. We will first treat the measurement 
on a phenomenological basis using the so-called positive 
operator- valued measure (POVM). In this formalism, the 
measurement is described by Kraus operators^ which 
we phenomenologically assume to have the form 

k[I] = [ 7:)07-e■''''*[*'^(*^(^"'^(*^"■^(*))/'="'^'(*)/^ . (28) 



Here T denotes time ordering and 



lR{t) ^ / dxi{x,t)g{x - xr) 



(29) 



is current operator averaged over a spatial region near 
the point xr > 0. By virtue of Naimark's theorem^S.^ the 
completeness relation / DIK''[I]K[I] = 1 is sufficient to 
guarantee that the measurement corresponds to a usual 
projective measurement in some extended Hilbert space 
including the detector. 

Additionally, we assume as usual that the detector and 
the system are initially uncorrelated.^ The spatial res- 
olution is taken into account by a convolution function 
g{x) = /2Aa; j ^/^t: ts.x ^ which is parametrized by the 
spatial sensitivity Aa;. To ensure, that the current on one 
side of the junction is measured, we also take Ax <C xr. 
The Kraus operator contains the parameter r, which 
plays the role of a measurement sensitivity. Varying r, 
the measurement changes from a weak, but nondemol- 
ishing measurement in the limit r — to a strong pro- 
jective measurement for t — > oo, which however yields 
not the expected result due to a large disturbing noise 
of the detector. For practical reasons, we convert dis- 
tances into times, namely, t„ — l^x/vn and t„ — xr/vu- 
The physical meaning of t„ is therefore the time, which 
an electron in the lead in channel n effectively interacts 
with the detector. On the other hand, i„ is related to 
the time-of-flight between the scatterer and the detector. 

The probability density of a given time trace for the 
current is given by a POVM^i 



p[/] = Trpift[/]if[/], 
The generating functional is defined as 



S[x] = In ^cxp 



(30) 



(31) 



I J dtx{t)I{t)/e 

where the average is defined as 

(...)^ = j DIp[I].... (32) 

Using the Kraus operators, one obtains for the generating 
functional 



S[x] j D<j) e 



S[x,</>]-/rft(20^(t)+x'(t)/2)/r 



(33) 



where is the standard Keldysh functional, defined 

as 



^(X(t)+20(t))4(i) 
ze 



S[xA] = lnTr|/5Texp 
X Texp 



^(X(i)-20(i))/fl(i) 



(34) 



The measure D(f> is scaled to keep 5[x = 0] =0. 

To calculate averages we need the transformations pre- 
sented in Appendix C. For future convenience we shall 
denote s±„ = s ± t„ and 



,{s) = e{s)h'^{s), K(s) 



727^7 



(35) 
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The mean current is 
SS 



{I{t))p = -ze 



Sx{t) 



(36) 



Hence, the conductance is robust against detector back- 
action as it does not depend on r. In the hmit of no-delay 
measurement (t„ — >■ 0), the measured current follows the 
(time-dependent) voltage 



(37) 



The correlation function is 



{SI{a)6I{b)), = 



e^P{a,h) (38) 



for SI = I — (I) p. We are interested in frequency scales 
of the current fluctuations in the regime w ^ t~^- Note, 
that the scale 1 /t„ plays a role of the maximal bandwidth 
of the detector. This means that the detector instanta- 
neously measures the current, which is a reasonable as- 
sumption also in typical experiments. 

In this case (see details in the appendix C) the noise 
measured at the detector contains several contributions 
and can be written as 

P{a,b) = -S{a-b) + Pr{a-b)+Po{a-b) + P^{a,b). (39) 

T 

The first term is just a Gaussian white noise of the detec- 
tor, which one would also expect classically. The Fourier 
transforms of the second term can be written as 



47r2 



^(2z0„(c^)+c^i?„e2'"*") 



9„(w)~-/i„(0)-iw/2 



(40) 



with 9n defined by ([55]) . Therefore, Pr has a frequency 
dispersion solely determined by the properties of the de- 
tector. It is independent of the bias voltage and can, thus, 
add an unknown contribution to the voltage independent 
background noise. Note, that it can be also made negli- 
gible by a suitable choice of the detector parameters, 
The third term has the form 



Po{u) - 



w{uj) 



(r2 + i?„(l-cos(2c^i„)). (41) 



Here w{uj) is defined by Eq. ([23l) . This term is also in- 
dependent of the bias voltage and can be neglected at 
sufRciently low temperatures. . 

Finally the most interesting contribution to the noise is 
the voltage-dependent part Pgi which we can write with 
the help of fell and Pe(a, b) = ^„ P„(a, b) as 

Pn{a,b) = RnTn-^^ — ^ COS <^(a_„, 6_„) (42) 
TTi n[a — bj 



with Tn{t) = exp (t(/i„(0) - /i„(i))/4) and 



ip{a,b) = $(a) - $(6) 



dteV{t)/h. (43) 



Still, this result can be applied to arbitrary fre- 
quency after multiplying Pn{uj,uj') by the damping fac- 
tor e~('^''+"'^)'^"/2. The noise (|42]) contains the usual 
symmetrized quantum noise in agreement with existing 
results The factor F results from the detector back- 
action on the current during the measurement on a time- 
scale set by the phenomenological parameter r. Note 
that T{t — >■ 0) — 1. A high voltage V leads to a strongly 
oscillating term cos((y9), so that P„ effectively probes the 
short time-scale t ^ h/eV. Hence, in the limit of high 
voltage we can take r{h/eV) — 1 and the noise is inde- 
pendent of the detector's backaction. To illustrate this 
effect more clearly let us take V = const and i„ = 0. In 
this case 

P„(w,w') = 2TrS{uj + uj')Pn{oj), 



w{t) 



cos{eVt/h). (44) 



>r„(0 

In Fourier space we obtain the convolution 

/d 
^P„T„F-i(a.-a) (45) 

X [w{a + eV/h) + w{a - eV/h)] . 

Due to the symmetry of T{uj) = F(— w) and the fact 
that r~^(w — ^ oo) vanishes, we can replace F~^(w) — > 
2ttS{uj) at high voltage and get the shot noise P„ = 
RnTn\eV\/Trh. Let us find the backaction corrections to 
the shot noise at zero frequency (a; = 0), 



/I 
^RnTMo')r-\a + eV/h). 



(46) 



For T <^ Tn, have 

r;;\a) = (27r - V^T/ATn)S{a) + Te'^'^^/A (47) 

and finally at zero temperature {w{a) — \a\) we get 

P„(0) = RnT4\eV\/Trh + qi\eV\Tr,/h)T/87:^T^] (48) 

with q{x) — e^^ — 2x dz . 

To illustrate the effect of the backaction due to the 
detection process we show in Fig. [1] the low- and the 
high-frequency noise for different detector parameters r. 
The upper panel shows that the effect of the detector is 
strongest for small voltages eV < ?i/t„. We can inter- 
pret this as follows: The additional noise added by the 
detection is dominant as long as the electrons flow with 
a rate smaller than the inverse interaction time with the 
detector < l/r„ through the contact. The influence of 
the backaction becomes negligible if the electrons flow at 
a higher rate, so that they do not feel the backaction and 
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FIG. 1: The voltage-dependent noise P„ at dc voltage and 
zero temperature for zero frequency cj = (a) and a finite 
frequency lj = f/r„ (b). In both cases the additional noise 
due to the backaction of the detector leads to drastic changes 
in the voltage dependence for eV < I/tu- The magnitude 
of the noise increase depends on the ratio t/t„. For larger 
voltages the usual shot noise is recovered. 



the full shot noise is recovered. A similar picture holds 
if the detector is sensitive to the finite-frequency current 
correlations. The signature of the quantum transport at 
huj = eV is gradually smeared out if the detector becomes 
slovifcr, viz. T > Tn. Interestingly, the noise first decreases 
vifith a finite voltage. Finally we note, that Fig. [T] shows 
numerically that the correction is small even at r ^ t„. 



A. Many detectors 

We can generalize the Kraus operator (|28| to the case 
of many detectors 

K[I] = [ D(t)Te^ '^'^^Aii'pA{t){iA{t)-iA{t))/e-4,\{t)/TA) _ 



(49) 

In particular, we take two independent detectors placed 
on the left and right hand side of the QPC. They are 



described similarly as in the single detector case, namely 



lA{t) = J dxI{x,t)gA{x - xa), for a = L,R. (50) 

Here we assumed xji > and x^ < and gA{x) — 

c-'^^/'^'^^^a/V^ttAxa with AxA < \xa\- The latter as- 
sumption means that the measuring regions of the two 
detectors do not overlap. We also denote the interaction 
time of the detector with the electrons as r^„ = Axa/vji- 
The conductance p7p remains unchanged but the noise 
can be now measured in four different ways. Namely, we 
can define the four correlators 



{6lA{a)6lB{b))p = - 



SxAia)SxBib) 



= e^PAB{a,h) 
(51) 

with A,B = R, L. 

We shall consider only frequency scales uj ^ Vn/\xA\, 
so that we can use t„ = in all derivations. Then, sim- 
ilarly as in the previous section (details in the appendix 
C), we obtain expressions for all correlators of the form 

PABia,b) = —SABS{(^-b) + PAB{a-b) 



TA 



+Po{a^b) + Pe{a,b). 



Here, Pq is defined by Eq. ([IT]) with t„ = 0. 
noise Pe is defined by Eq. but with 



r„(t) = exp 



M-hAn{t))/4. 



(52) 
The excess 

(53) 



Note that the excess contribution is the same for all cor- 
relators and depends on the parameters of both detec- 
tors. The voltage independent contributions differ for 
the auto-correlators and the cross-correlators. We have 



Paa{uj) = PrA{uj) + 



47r2 



(54) 



where P^a is defined by (|40p with t = ta and t„ = 
and B = L,R for A = R,L, respectively. Finally, for 
A^ B wc find 

Pab(-)-£jE^"E (55) 

n m 

[TA{2i9A7n{uj) - i?mW) + Tb {-2i9*g^{uj) - i?,„w)] . 

The possibility to measure several independent correla- 
tors has interesting consequences. In the case of a single 
detector one could in practice measure only the voltage- 
dependent contribution of the noise as the offset noise 
has generally an unknown value and is subtracted. How- 
ever, the use of two independent detectors helps to es- 
timate the background noise in the auto-correlation sig- 
nal. Hence, comparing cross- and auto-correlations we 
also can get a rough estimate of the offset noise. Fur- 
thermore, in the limit tl,r — we get Plr = Po + -Pe- 
As Pe is independently known from the auto-correlation 
measurement, we can get information about the voltage- 
independent part of the noise Pq. 
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B. Higher cumulants 



rections to (l57l) of the form^— 



Now we consider the effect of our detection scheme 
on the third and fourth cumulants of the current fluc- 
tuations. While these are harder to measure than noise 
correlations, they contain a tremendous deeper informa- 
tion that the current. It should be noted, that the third 
cumulant already at zero frequency contains a non-trivial 
ordering of the current operators^ and we can expect a 
similar non-trivial effect of the detection scheme. Quite 
generally the third and forth cumulants are defined as 

{{ABC)) = {dAdB6C), (56) 
{{ABCD)) = {SASBSCSD) - {SASB){SCSD) 

-{SASC){6B6D) - {SASD){6CSB) . 

Here SX — X ~ {X) are fluctuations of some observ- 
able. In the following these will be current fluctuations 
operators at different times. 

We now apply out model of a detector, which is 
parametrized by a Kraus operator to calculate the cu- 
mulants. Unfortunately, the general expressions are too 
cumbersome to be shown here and we discuss only a limit- 
ing case below. Hence, we take the limit |TnU;|, t/t„ ^ 1, 
where w describes relevant frequency scale of the mea- 
surement. In this case the only effect of the detector is a 
large white Gaussian offset noise l/r. It only adds as a 
constant to the second cumulant, while higher cumulants 
are unaffected. 

The derivation of the third cumulant is given in the 
appendix D. The resuh for C3(a,6,c) = ((/(a)/(6)/(c))) 



C3 (a, 6, c) = ^ I RnTn {Rn - T„) X 
^^^°r"^ (5(a - b)5{a - c) + V f ds P^{a, b; s) x 



a{abc 



[(2T„ - l)w{s - c_„) + w{s - c+„)] tt/2 



(57) 



where the summation is over all permutations of the set 
abc. At zero flight time and either zero temperature and 
frequency or tunneling limit only the first term survives. 
The last part contains the so-called noise susceptibility-'^* 



P„(a,6;s) = 



hdPn{a,b)[V] 
edV{s) 



(58) 



generalized here to an arbitrary time-dependence of the 
bias. The result (|57|) in the limit of zero flight time has a 
somewhat simpler form than the existing results The 
identification of the noise susceptibility in the third cu- 
mulant has important consequences on its detection. A 
general detector, which is equivalent to an electromag- 
netic environment, and weakly coupled here, gives cor- 



C3"(wi,W2,cJ3) = X! y f^wP(wj,Wj;'^) X (59) 

<j{ijk) 

[boiuj)w{uj)Siu} ~ ujk) + bi{ijj)P{-uj,uJk)/e^]/A. 

Here P{a,f3;s) — SP{a, (3) /SV{s) is the noise suscepti- 
bility. The function 69 represents the influence of envi- 
ronmental noise and 61 is the voltage-dependent feedback 
of the environment due to the noise of the system. They 
can in principle be modeled by an effective electric cir- 
cuit. In practice, the precise environmental circuit is not 
known, and the functions are determined by fltting some 
model environment.^^ 

An important case in practice is a constant bias volt- 
age, V = const. The expressions for the noise suscepti- 
bility and the third cumulant are simplified considerably. 
The noise susceptibility turns out to be 

Pn{a,f3;u) = R,,TJ{a + /3 + ij)e-'^*" X (60) 
[w(7 - eV/h) - w(7 + eV/h)]/u] . 

7— a,/3 

The noise susceptibility is shown in Fig. [5] for different 
temperatures. In all cases, it preserves the symmetry 
a o ±/3. At zero temperature, a non-analyticity occurs 
along lines \eV/h\ — \a\,\/3\, which is smeared out at 
increasing temperatures. 

The third cumulant also becomes much simpler. Due 
to the relation C3{a, f3,uj) — 2Tr6{a + /3 + uj)C3{a, (3), 
it effectively depends only on the frequencies a and /3. 
The delta function imposes the constraint a + P + uj — 
so only two of the three frequencies are independent. 
The cumulant is hence symmetric under changes a, /3 
uj — — f3. In the following expressions, we either use 
two independent arguments {a,/3) or three constrained 
{a,l3,uj). For a neghgible flight time |a|, |/3|, |a;| ^ we 
have 

C3{a,f3) = ^i?„r„(l-2T„F(a,/3))2e4T///i, 



F{a,l3) = 1- ^ u{uj)h/eV X 

[w{a - eV/h) - w{/3 + eV/h)]/4, , 



(61) 



where we sum over permutations of the constrained set 
al3itj . The temperature-dependence is encoded in a single 
universal function F, which does not depend on the chan- 
nel transparency. In agreement with previous results^li^ 
we find at zero temperature 

F{a,l3) = l-min{max{|a|,|/3|,|a-)-;3|},|ey/;i|};i/|ey|. 

(62) 

The zero-frequency limit is of course^ 

sinh U — U , , , 

^(»'")-^- V(cosh^^l) '^--^/^-^- ^''^ 



The frequency dependence of the function F is plotted in 
Fig. [31 The form of the plot is motivated by the symme- 
try of arguments and the constraint a + (3 + lu = 0. How- 
ever, experimentally a verification of the full frequency 
dependence is a desirable challenge. 
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FIG. 2: (Color online) The susceptibility P„(a,/3;cj) at dc 
voltage defined in (|6Up for different temperatures. The cross 
symmetry is preserved for all temperatures. 

The effect of flight times on the second cumulant in 
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FIG. 3: (Color online) The factor F{a, /3) defined in ((6T| for 
different temperatures. The hexagonal symmetry is preserved 
for all temperatures. 



Eq. (|4ip was only an additional offset noise, independent 
of bias voltage. In the case of the third cumulant the 
terms, which depend on the flight time, depend on volt- 
age, see Eq. ([57| . For :s> we get 



C3iuj,uj') = Fiuj,uj')C3iO,0) 



(64) 



so the cumulant drops to zero at large frequencies, con- 
trary to the zero flight time case, where it stays nonzero 
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and proportional to eV. The recent experiments^ does 
not show any reduction of the third cumulant at high fre- 
quency, so the frequency scale defined by the flight time 
must still be beyond experimental reach. However, there 
is an environmental correction of a similar magnitude. 
We will continue the discussion of the third cumulant in 
Sec. IV. 

Lastly, the fourth cumulant can be evaluated in the 
similar way. Here we present only the result at eV = 0, 



((/(u;i)/(c^2)/(c^3)/(^4))) = 

[w{u]l) + 'W{UJ2) — + Wa) — w{ljJ2 + W3)] 



(65) 

^ u{uj^)u{uj4) X 



where we perform summation over all permutations of 
the set cji, . . . , aj4. Due to hyperbolic identities, for tn = 
it reduces to 



E 



p^E T 



(66) 



which agrees with existing results at small reflection or 
transmission, s* One can in principle use this method 
to find analytic expressions for cumulants of any order, 
which is beyond the scope of this work. 



C. Tunneling and transparent limit 

In some limits the generating function S[x\ can be eval- 
uated exactly. For r — 0, we can make the separation 

S[x] = S[xA- I dtxHt)/2T, 

where S[x,(j)] is defined by (jC4[) . We are particularly 
interested in the tunneling limit T„ <C 1 and the trans- 
parent limit Rn ^ 1 up to terms of the first order in r„ 
(i?„).S'' With help of appendix E, we get for i?„ ^ 1 



S[x,0]^2NSo[x] 



ieV{t) 



n 71 

(67) 



and for r„ ^ 1 (tunneling limit). 

Six, 0] = ^ 45o [xTn] + ^Tnix] • (68) 



is the generating functional of Gaussian vacuum fluctua- 
tions and 



SAnix] 



ieVU) , . , 

CA Sm X[t+n jdt + 

Trn. 



dtdt' ""^^ ^^ A„ <j eAXAn{t) sinx(i'+«) (70) 



. X{t+n) X{t'+n) r^{t' ,t)+D{t' ,t)\xAr.] 

2 2 



The arguments depend on the transmission (reflection) 
probabilities and are defined by 2xAn{t) — xi^-n) + 
CAxii+n), ^R/T = +/^- The kernel in the exponent is 
defined as 



(71) 



D{a,b)[^]^ J ds J dtw{s-t)^{t) 



Re 



, , n(a — t) 



2tt 



coth + le 



h/keT 



One can see that the most significant effect of the flight 
times in the tunneling limit is Gaussian noise S'o that is 
growing with distance from the contact. This is because 
XT{t) — x{t-) ~ x{t+) vanishes at zero flight time but 
for large flight times the two parts will become indepen- 
dent and not cancel each other. The noise saturates to 
the same equilibrium value as in the transparent limit. 
We stress again, however, that there is no experimen- 
tal evidence of reaching that timescale so far. The non- 
Gaussian part Sau remains small as it is proportional to 

An- 

In the limit of a vanishing flight time i„ = the above 
formula for the tunneling case simplifies to 



dtdt 



idt sinx(i) eV{t)/Trn (72) 
,2w{t~t') . x{t) . X(t') 



■ cos ip{t, t') 



with w and ip defined by (|23|) and (|43l) . respectively. 

It is tempting to interpret the last result in terms of 
a counting statistics. Namely, we might identify terms 
g±jx(t) and 6^*-*^^*^/^^*-*^^* with a quasi-charge transfer 
of ±e at t and (±e/2,±e/2) at {t,t'), respectively. The 
fact, that in this expression half-integer charges appear 
has probably a similar origin as the half-integer charge, 
which appears in resonant tunneling. Hence, inter- 
preted as a quasi- Poissonian distribution we may identify 
"probabilities" according to 



Here, N denotes the number of modes 

.wit-t') 



5*0 [C] = - / dtdt 



mm, (69) 



S ~ ^p,(t)e^^xW 



(7=± 



Jcrxit)/2+ia'xit')/2 



(73) 
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By comparing Eq. (l73l) to Eq. (l72l) we can read off the 
transfer "probability" per unit time. These take the form 

p,{t) = aJ2TneV{t)dt/2nh, (74) 
Paa'{t,t') = aa'J^Tndtdt' '^''^^^^ cos (p{t,t'). 

n 

Unfortunately, these rates can be negative, so they can- 
not be interpreted as a probability. Only a combina- 
tion of and Pa-a-' would make sense, which actually 
happens in usual full counting statistics - valid at long 
times. ^ Hence, the generating function (17^ itself does 
not correspond to a measurable probability. Only after 
convolution with the Gaussian detection noise we get a 
real probability. 

IV. QUANTUM TAPE 

In section III we have introduced a measurement pro- 
tocol by means of Kraus operators. Here we go further 
and want to find a quantum detector that is coupled to 
a source and registers the time dependence of the source. 
Our aim is make a quantum tape that translates time 
dependence into spatial dependence. The tape interacts 
with the system at a fixed point in space. Then the tape 
moves far from the point of interaction with the source 
and afterwards the projection (reading) is applied. 

Our quantum tape will be a linear quantum wire or 
equivalently a massless Josephson transmission line. In 
this description we have the joint Hamiltonian of the de- 
tector and the system 

H = Hd + Hi + Ho- (75) 

Here 

H'^ = ^j ds{Ql{s)+il{s)), 

Ho^J2 f dx{v{t)ei~x)QrAx) (76) 

+ihVn[%l;\^(x)dx^Ln{x) - il'^j^^{x)dx^Rn{x)], 

+qn5{x)['4j\^{x)'iijBn{-x) + (x)'0L,-i (-o:)] | , 
Hi^-'^ j ds\{s)Qd{s)Q 

with A(s) = A(-s), Q ^J2nl dxfn{x/Vn)J2nQn{x), 

Qn = X! ^'^\n{^)^An{x). (77) 
A=R,L 

The bosonic operators satisfy (see also the fermionic rep- 
resentation in appendix F) 

[h{s),h{s')] = [Qd{s),Qd{s')]^Q 

[/d(s), Qd{s')] = -le^dsis - s')/n . (78) 




FIG. 4: The system (quantum point contact) coupled to a 
one-dimensional wire (a) and possible equivalent electric cir- 
cuit (b). Ammeters, switched on at instant time, measure 
spatial profile of the current. 

We used here time units for the Hilbert space of the de- 
tector (Fermi velocity Vd = 1). The functions A and /„ 
define the coupling between the detector and the system. 
Their arguments are in time units and should be nonzero 
only on one side of the junction. In principle it should 
be the total charge, i.e. /n(s) = 9{s), but we have to al- 
low finite range of the coupling. The setup is depicted in 

Fig. m 

We also define Heisenberg operators A{t) — 
U^t)Au{t) with U{t) = Te-fo mt')dt'/zh^ rpj^g — i-ij^i 
sity matrix is p — pdPs with pa oc e^^'^f^sTa a = d,0. 
The current operator is defined / = e.ip^f^^'Ln ~- L ^ 
R. We define the auxiliary current operator 

Iis) = iidis)-idi~s))/2. (79) 

Now, performing a spatially resolved measurement of 

I{s,t) at s ^ (beyond the support of A(s)) we get 
the probability density functional 

p[i]^TTp5[i - i]. (80) 

The generating functional S[x\ = ln(e^'-^'^^^"'^*^''^*'''^'*/'^) in 
the limit of weak coupling (neglecting terms ©[xA^]) will 
be given by (see appendix F) 

S[x] = S'[-x*A,ix*A*u;d/2]-y da;|x(a;)pWd(a;)/167r^ . 

(81) 

The convolution in the arguments is defined as (a * 
A)(w) — a{Lu)X{uj), x{uj) — ~iujx(y>)- and 

5[x,0] =liiTrpore/'"'*W*)+2*(*))<3«W/2^ x 

J- gJidt{xit)-24>it))Qoit)/2e^ ^g2) 

Here we used the interaction picture, in which Aolt) — 
U^{t)AUo(t) with Uoit) = refo^°^*')'^*'/^^. We get the 
expected detection picture - the measured signal I{uj) 
is proportional to A(a;)/o(cLi) with Io{t) = dQ(t)/dt = 
i[Ho{t),Qo{t)]/h. The contribution from Iq is much 
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weaker (A ^ 1) than the internal Gaussian noise of the 
detector - the last term in (|5T|). To evaluate ([5^ . we 
can apply the FCS formalism developed in Sec. Ill and 
App. A for Hj = J dxJ2nQfii^Wsn{x,t) with the bias 

Vsn{x,t)+V{t)0{-X), 



Vsnix,t) = h(t){t)fn{x/Vn)/e. 



(83) 



The averages measured at the detector are a combi- 
nation of its own noise, the noise of the system and the 
response of the system due to the detector. The average 



(84) 



(U 
DC 

CM 

3; 

CM 
05 

10 

05 

DC 




2.5 3 



with 5„ = 5„(1 + 0[g'^]), 5„(a;) = -ia;A(w)/„ja;). 
The observed fluctuation spectral density of / is 
C2(w,w') = {5i{<jj)5i{ijj')) = 27re2,5(w + a;')P(w) + 
e^5„((u;)5„((u;')P„(a;, w'). It contains a voltage inde- 
pendent contribution 

P[uj)=Wd{uj)lA^-Y,RnTnW[uj)\gn[L0)\^lTT (85) 

n 

+ ^[i(;(w) - u;rf(w)] [|5„(w)P - i?„Re5,'.(w)] A , 

n 

where it;(w) = a; coth(fta;/2fcBTo). The second part is 



, I. w(t-t') 
Pn{t,t') = fl„r„ \ ' cos 

2sin^(wt) 



rf,£ZM, (86) 



lnr„(t) ~ duj 



(27ra;)2 



:i«rf(a;)|5„(c.)|^ + 0[ff4]]. 



The excess noise - i.e. the voltage-dependent part of C2 
- is contained in (|86p . Note that this result is valid for 
an arbitrary time-dependent bias voltage. In the limit 
of weak coupling, viz. gn <C 1, we can take T{t) ~ 1 
and the excess noise agrees with the symmetrized quan- 
tum noise. It is a tiny contribution on top of a huge 
background noise P(w). The background noise is also 
seen in experiments but it cannot be distinguished from 
the amplifier's noise. Measurements of the voltage de- 
pendence reveal only the excess noiseJ^^ In the limit 
of strong coupling, (7„ ^ 1, the noise is heavily affected 
by the detector's backaction. The excess noise is affected 
through T(t) > 1 and gets a significant correction, sim- 
ilar to that in Sec. Ill, which vanishes for high voltage 
since T{t ^ 0) = 1. 

The situation is different for the third cumulant. In the 
limit of weak coupling we keep only terms to the lowest 
order in gn, which gives 



C3(a;i,a;2,W3) = {Si{uji)Si{uj2)Siiuj3)) 



(87) 



^eV{u}i +UJ2 + W3) 



+ ^ (2T„ -f g„(wfe) - l)P„(wj,a;j;a;fe)u)(wfe)/4 

<T{ijk) 



FIG. 5: The high-frequency third cumulant of current fluctu- 
ations in the tunneling limit at zero temperature for different 
values of g = g„. The average current is Iq = 2e^ TnV/h, 
we took \Q\ <g {uj,eV/h) and assumed a mode- independent 
coupling g = gn- Note that Reg < at = and g = at 
ksTa = 0. 



for the susceptibility P„(a, /?; s) defined by ((58)) . (7„(w) — 
ei^^i'^)[l-Wd{uj)/w{uj)] and 7?„(cj) = -2arg5„(a;). The 
inner summation runs over all permutations of the set 
123. Note that g„ leads to corrections that are indistin- 
guishable from the infiuence of an environment described 
by ([5^ (up to factors g„). The detection dependent part 
g„(cj)((7„(a;) — l)/h from ([57]) can be effectively absorbed 
into bo in ([5^ . Hence, the detector can be interpreted as 
an example of some environment. 

Another interesting observation is that all cumulants, 
except the first term in (1851) . vanish at equilibrium be- 
tween detector and the system, namely for eV — and 
ksTo = ksTd- On the other hand it is obvious because 
no information transfer is possible at equilibrium as the 
entropy is already maximized. An eflacient detector can- 
not be in equilibrium with the measured system. 

Let us finally consider two limiting cases. Analogously 
to (HH) and for fcsTo = and Q < (w, \eV\/h) we 
get from dHZ]) 



C3in,u;) = ^5„(^])|.g„(c^)p2i?„^„ — 



l - 2r„ -f (2r„ - 1 + Req„(a;)) min 1 



h\uj\ 
W\ 



This result is plotted for different g„ in Fig. [5] in the 
tunneling limit, r„ ^ 1. We see, that g„ induces a step in 
C3 for eV = huj. For arbitrary ksTo and \uJk \ ^ \eV\/h, 
we have 



C3(W1,W2) = 2RnTn ]"[ gn(^fc) X 

n k 

1 - 2T„ + ^2T„ - 1 + ^ 9n(^fc)/3^ (1 - F(0, 0)) 

with F(0,0) given by in agreement with the envi- 
ronmental correctioUiiS We see, that the behavior of the 



(89) 
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third cumulant is different from the noise. Even in the 
weak coupUng hmit it may get a significant quantum cor- 
rection due to the detector. 



A. Many tapes 



then classical while stih gaining quantum information. 
Such a situation is close to the POVM model from Sec. 
Ill in the limit of weak coupling because the detector's 
noise (last term in ([Ml)) becomes white and decouples 
from the system's signal. 



We have seen in Eq. (|8T|) that the outcome of the mea- 
surement is always correlated with fluctuations in the 
tape as the second argument of S is non-zero. This can 
be avoided by means of many tapes, namely 



(90) 



with 



[Lis), Qb{s')] = ~ie^5abds5{s - s')/n, (91) 

[ia{3)Jb{s')] = [Qa{s),Qb{s')]^Q. 



We also redefine 

/(.'^)-E(^^(^)-^^(--'^))- 



(92) 



The resulting generating function (jFll|) has the new form 
in the limit of weak coupling, 



S[x] = s 



-^X* ^c,^iX* ^c* Wc/2 



dw|x(t^)/47rp^Wc(a;). 



(93) 



The results for the cumulants differ from the single tape 
case as follows. In the definition of gn below (|84|) we have 
to replace A J2c ^c- In (|85l) we change the first term 
Wd Z^c^c- We also change Wd Yl^c^cWc/ J2c-^c 
in the second line of (|85p and the definition of qn below 
(|87p . This change makes in principle arbitrarily large q„ 
(positive or negative) possible even at = 0. On the 
other hand a simple model of environment'^ - can give ef- 
fectively (?„ < 1 when comparing ([87l) and (|59)) . while the 
case (?„ > 1 has been observed in the recent experiment 
The latter suggests that the actual detection scheme or 
model of environment is more complicated. 

Let us consider a very special case of tapes at high tem- 
peratures 1/AcT ^ ksTc/h 3> 1/t, where r is a certain 
time resolution so that the measurable range of frequen- 
cies is |™| < 1. If J2c ^cksTc = (but Y.C i- 0) then 
we can write 

S\^=~S\-Y,X^\,A- f dL.|x(^)PE«;c(^)/16^^ 

c c 

(94) 

so that the system-detector correlation cancels (0 in the 
second argument of S). The detector model becomes 



V. CONCLUSIONS 

We have presented detection models of time-resolved 
quantum detection of current in mesoscopic junction. 
One can make the measurement by means of either a 
partial projection (Kraus operator ~ POVM) or a full 
projection on a weakly coupled transmission line - quan- 
tum tape. It is difficult to separate the backaction of the 
detector from the signal of the system. Nevertheless, the 
generality of the presented method makes it applicable 
to a wide range of types of measurement. 

Applied to a quantum point contact, both models give 
the expected expressions for current and noise if the back- 
ground noise is subtracted. The voltage-dependent part 
of the noise is independent of the detector in the weak 
coupling limit. On the other hand, the third cumu- 
lant may contain significant voltage-dependent correc- 
tions even in the weak coupling limit. The correction 
is indistinguishable from the effect of the environment. 
Hence, the high frequency measurement in mesoscopic 
junctions always contains a detection-dependent contri- 
bution. Experimentally, an independent determination 
of the coupling parameters would help to distinguish be- 
tween the effects of the detection process and the envi- 
ronmental backaction. 

We have also derived an analytical tool to calculate 
frequency-dependent cumulants for a mesoscopic junc- 
tion with energy-independent transmission. An impor- 
tant result is to identify the noise susceptibility in the 
expression for the third cumulant at high frequency. The 
method works well both in time and frequency domain 
and can be extended to other complicated correlations. 
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Appendix A 

In our model of the quantum point contact, we make 
the approximation of a linear spectrum and a short range 
scattering potential. Both assumptions lead to singular 
behaviour and we have to regularize some terms in the 
Hamiltonian and describe the details here. 

The short range potential scattering is introduced in 
Eq. ([5]) and we would like to relate this term to the 
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scattering matrix. One of the possible regularizations of 
® is 



(Al) 

with the function Zn{x) — >■ qnO{x). For |fc„a;| ^ 1, the 
wave eigenfunctions satisfy 



ihdx1pLn{x) + z'n{x)ll>Rn{x) = 0, 
-ihdxtpRn{x) + z'^{x)'lj}Ln{x) = 0. 

We can change variables x ^ z{x), 



(A2) 



(A3) 



which gives the solution 



c+e 



+ C-e 



'ic+e"'" " + ic^e 
By taking ipLnizn = g™) = 0,ipm{zn = 0) = 1, we get 

1 



(A4) 



cosh{qn/h)' 
fn = ^Ln{zn = 0) = -i tanh(g„ . (A5) 

The phase factor in i„ is completely irrelevant in our 
approximation. However, it can be added by a constant, 
symmetric bias potential. 

For the future convenience we shall use a little different 
regularization of the second term in ([5]) , namely 

^ / dxz'^{xlVn) '!i^''^^{x)'^Rn{-x) + '!i}''j^^{-x)'^Ln{x) 
n •' 

(A6) 

It gives tn = cos(q„/fi.) and f„ = i sin((7„/fi,). The regu- 
larized Hamiltonian (fT2l) takes the form 



H = Ho + Y, [ 2z'„is)hn{s)/e 

+ ^ y" V{±SVn, mon{s) ± /ifi(s)] (A7) 



when acting on the density matrix standing to the right 
of these operators. Next we construct the operators 



dk 



?/^L(a;)V^Lfi(a;) = ^ 



-ikx 



A{x) 
A{k) 

The last operator needs to be regularized as 
i(fc) 



i(fc). 



^'^ --9{q)i'lni<l)^Lniq + k), 



(AlO) 



(All) 



where g{q) — 1 for \q\ <C A and g{q) — ^ for g A, 
with A denoting some cutoff larger than all relevant wave 
scales. Then 



Hence, 



[i(fc),i(fc')] - / ^i^l^{q + k + k')i^Ln{q) X 



(A12) 



5(9)(5('Z + fc')-5('Z + fc))- 



(A13) 



For |q|, |fc|, |fc'| ^ A the above expression gives zero be- 
cause g{q -f- k') — g{q + k). Also for \q\ >• A, it zero 
due to q{q) — 0. For q ^ —A we get again zero due to 
(|A9|) . The only nonzero terms come from g ~ A. We 
can use here also (jA9|) but only after changing the order 
of the operators "0 and "0^- The net contribution is the 
anticommutator for q > 0. Hence, 



[i(fc),i(fc')] = 

^5{k + k')g{q){g[q + k')-g{q + k)) 



(A14) 







27r 



dq 



= <5(fc + fc') / ^g{q){g{q - fc) - .9(g + k)) 
Jo 2tt 

S(k + k') [ ^g{q)g{q - fc) = k5{k + k')/2n. 

In the last line we used the fact that g{q) — 1 for q,k 
A. By taking the Fourier transform we get 



[A{x),A{x')] = idj{x - x')/27r. 



(A15) 



To derive ([T^ and P3|) . we recall the standard bosoniza- 
tion scheme for a ID system of noninteracting fermionsj^ 
Let us define 



V'Lfi(fc) = 



dx 



e'^^^Lnix). 



(A8) 



We assume that the excitations of the Fermi sea are finite 
in the sense that levels deep below and high above the 
Fermi level are always occupied and empty, respectively. 
This leads to 



iJ'Lnik -OO), 0lft(fc ^ oo) ^ 0, 



(A9) 



Appendix B 

The current operators defined in (|lll) can be found 
exactly by a solution of their respective equations of mo- 
tion. We have the following Heisenberg equations ([M)) 
for (dB), 

Dion{s,t) = e^d,V+{sVn,t)/ATT, (Bl) 
Dh^{s,t) = -2z'„{s)i2n{s,t) + e^d,V^{sVn,t)/4TT, 
Dl2n(s,t) = 2z^j(s)/ift(s,i) - eV-{sVn,t)hn{s,t), 
Dhn(s,t) = eV-{sVn,t)i2n{s,t), 
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where D = h{dt - ds) andV±{x,t) = V{x,t) ±V{-x,t). 
For s — >■ +00 we have asymptotic equilibrium operators 



for j — 0, 1, 2, 3 and ts =t + s. The operator Iq decouples 
from the other currents and has the general solution 



4ft(s,i) = /ofi(is) + / ds'V+{s'v„,t + s-s') 



(B3) 



For s > we have 

hn{s,t) = im{ts) + -^a,$„(s,t), 

47r 

hn{s,t) ^ /2fi(ts)cos$„(s,i) - /3ft(ts)sin$„(s,t), 
hn{s,t) = /3^(ts)cos$„(s,i) +/2n(is)sin$„(s,t), 
e /"^ 

^n{s,t)^-j ds'V-{s'Vr,,t + S-s') (B4) 

and for s < 

knis.t) = i^,-^{t s) COS <^>+n{s,t) - /^j,(t,)sin$+„(s,i), 
l3n{s,t) = /^,-,(ts)cos$+„(s,t) +/°-^(t,)sin$+„(s,t), 
$+„(s, t) = $„(s, t) - $„(i) = $„(0, i), (B5) 

with the boundary conditions 



ilnis) = [Tn - i?„)/lfi(0+, S) - 2v'i?„T„/2n(0+,.s), 
= (T'n - i?„)/2fi(0+, S) + 2 /R^/ir (0+ , s) , 
4(s)=/3(0+,s). 



In particular for s < 

/lfi(s,0 = {Tn - Rn)hn{ts) 



(B6) 



+— 9,($„(s,t) - 2i?„$„(t,)) - 2v/i?«r„J*(t,) 
47r 



with 

i|(s) =l2fi(s)cos$„(s)-/3ft(s)sin$„(s). (B7) 

The great advantage of the above equations is that the 
Heisenberg current operators can be represented by com- 
binations of equilibrium operators. The same applies to 
averages and facilitates the calculation of higher cumu- 
lants. 



Appendix C 

To evaluate the noise correction due to the detetctor 
in the POVM model we employ the the formalism devel- 
oped previously. We have to perform a number of path 



integrals to find the final analytical expressions for the 
noise. We start by splitting Eq. into two parts 



(B2) lR{t)=T.[^m{t)+Imit) 



(CI) 



Imi^)^ J dxg{x + XR) Im{xn,t) - Iin{xn,t) 
Imit) = I dxg{x-XR) -ion{xn,t) - hn{Xn,t) 



with Xn ^ x/vn- To simplify (j34p . we will use the com- 
mutators 



[/fl\W,/r(O]=0fort>t' 



(C2) 



ZTT 

with hn defined by (I55|) . Using the Baker-HausdorfF for- 
mula e-^e^ = gA+Sg[A,s]/2 .^g gj^^ 

idt 
2^ 

-idtdt' 



re/^«(*)^«(*) = J]exp 



exp / ^^(t)4"nWexp 



Using similar transformations we get 



mimit) X (C3) 

mat'Wnit-t'). 



(C4) 



where 



C/i" = exp 
exp 



rrout 



idt 
2^ 
idt 



e 

'out 



f7^f, = exp / 



dtdt' 



X{tm'Wn{t - t') 



(C5) 



with 6n defined by (1551) . Next, to find the noise (P^ . we 
need the help of the relation 



ZJ^exp 



202 (i) 



(C6) 



which gives 

P{a,b) ^ 5{a-b)/T + 

I dtY,^()nia~t)eUb-t)+ I i^0e-^^^'W 



E 



2di 



:<^(i)e„(a-<)Tr/54^(6) 



o6 



(C7) 
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and SA = A — TipA. Finally, we get 



Mt) - + / V2dt'(j){t^n + t')K{V2t') (C8) 



where J is defined by (jB7l) with $„ replaced by </>„. 

The case of many detectors (Eqs. (|49)) - (|55p ) is analo- 
gous, 



iA{t) = Y.[rXn{t) + iTnit) 



(C9) 



I°Anit) = j d^ai^ + \xA\)[<^Al0niXn.t) - im{Xn,t) 

iAnit) = / dxg{x - |a;^|) -eAlonixn,t) - im{Xn,t) 



with eA — sgn xa ■ In Eq. (jC4p we put 



idt 



idt 



idt 



U^n = exp / — ^ 



/ '^Y.^AmA{t')eAn{t-t'). (CIO) 



We finally obtain 



PAs(a, 6) = 5{a - b)5AB/TA + 



D(l)e 



+ E / ^Mt)OAnia - t)TTpit^ib) 

nm 



(Cll) 



and 



eAlonit) - (T„ - i?„)/i^(t) + 2^R^Jf^{t), 

(t>n{t) = ^{t) + ( V2dt'^dpA{t + t')hAn{V2t'), 

A 

with B = L, R ioi A = R, L, respectively. 

Appendix D 

The third cumulant ([57)) has the form 

((/(a)/(6)/(c))) = ^TrpF4a,&,c) . (Dl) 

n 

The operator in the last equation has the form 

"1 



Fn{a,b,c) = ^ 

aiabc) 



i(/i„(a+„)/-(6_„)/-(c_„) 



+/*(a_„)/?(6_„)/ifi(c+„) 



(D2) 



Here a summation over all permutations of the set abc is 
assumed and 



lt{t) = {Rn - T„)Im{t) + 2^R^TJt{t). (D3) 

In these equations J is defined by (1^ and $ by (P7|) . 
Using the commutation rules (IT4|) and (|2T|) we gelj^iiS^ 

((/(a)/(6)/(c))) = ^ i?„T„e3 ^ sin (/9(a_„, 6_„) 

n CT(a6c) 

) +(7(a_„,c+„,6_„)]/2 + 



(i?„ - Tn)^^^^—-^S{a - b)S{a - c) 



(D4) 



The Fourier transform of q{a, c, b) has the form 

qia, -/,l3)^iS{a + p + 7)u(7)(w(a) - w{l3j), (D5) 

where u and w are defined by (j23p . The function is 
defined by Eq. ((43|) . One can also find that^^ 



q{a, c,b) ~ w{b — a) / ds w{s — c) 

Jb 

sinhfB - A) 
Re ^-^ — X 



sinh^(B- A + ie) 



Re- 



(fcsT/ft)3 



2 sinh(C - A - ie) sinh(S - C + ie) 

for X = irkBTx/h, x = a,b,c and e — )■ 0. Special care 
must be taken at a = b, b = c, c = a. Note, that the 
function g(a, c, 6) is not cyclic because 

q{b, a, c) - (7(a, c, 6) = 96[(5(a - 6)(5(c - 6)]/27r. (D6) 
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Appendix E 

In the limiting cases <C 1 or T„ ^ 1), it is possi- 
ble to find expressions for the generating functional (I7ni) . 
Here we present some details on the derivation. The fol- 
lowing operator expansion formula is useful, 



1 i-l-x 
"'0 



(El) 



We use it for A — [Rn — r„) / dt ix{t+n)Iin{t) / e and 



B — 2\/ i?„T„ dtix{t+n)J^{t)/e. (E2) 



What we need is the term of the second power in B and 
an algebraic identity 



1 /-i-j/ 

dx / dy ^7jexp[i{x — 1/2)^ + i{y ~ l/2)ri 

"'0 

2 sm - sm — z sm — i- sm — - — . 

2 2 2 ^-rj 2 



(E3) 



We also need several auxiliary operator identities. Let us 
define 



and 



E^[^] = exp - / idt at)hnit)/e (E4) 



5oK]=hiTrp^|[e] 



(E5) 



for an arbitrary function ^. One can show using (jl4p and 
(PT|) that 



dt 



= iu{u)TTp[El[^],h^{iu)] (E6) 



and 



So 



[EmJmico)] = EM I (E7) 



,t6So[^] _ w{uj) 



dt e*"*C(0 (E8) 



sm 27r 
and we finally obtain (j69p . Another useful property is 

I±n{t)En[^] = Enm±n{t)(^^'^^'^ . (E9) 
We define the operator 

Dn{a,bm = Enm+n{a)i-n{b)En[^] (ElO) 



with i±n{t) = i2n{t) ± ii3n{t)- We Can show that 
,STvpDn{a,bm 



Idt e 
and 



sm 



= iu{oj)TTp[Dn{a,b)[£,]Jin{uj)] 

(Ell) 



[D^(a,6)K],/ifi(c.)] = A.(a,6)K] 



2tt 



'm 



(E12) 



This gives 



^^ TrpA.(a,6)[g] ^ 
Trpi?fi(a,6)[0] 
where D was defined in Eq. (ITT]) . Finally, we get 

^[X, 0] = ^ 2(5o[xo„] + ^o[xi«]) + 2 ^ 5„[x], (E14) 

n 71 

where 2xon(i) = x{'t-n)-x{t+n) and 2xi„(t) = x(*-n) + 
(T„ — i?„)x(i+n)- The lowest order correction is 



^n[x] ^ 
Xi.t+n) 



dt{R,,TJ4^>^ 



5Xln{t) 

sin((r„ Rn)xit+n)) 



ieV{t) 

2TTh 



Xi^+n) ~ RnTn 



Sin((r„ - Rn)x{t+n)) 



(r„ - i?„)2 



^.^(T„-i^„)x(t,„)^.^(T„-i^„)x(^VJ^ 

2 2 ^ ^ 

We stress that the above formula are exact up to first 
order in T„ ^ 1 or i?„ ^ 1. After taking the respective 
limit higher order contributions should be disregarded as 
we would need also B^-terms. 

Note also that Ion = {Rn - Tn)hn + 2Vi?«r„/2fi sat- 
isfies the same commutation rules with itself and the 
Hamiltonian as Im. This leads to the identity 



= J dtdt' w{t - t') X (E16) 
X(0sinx(i)-4sin^sin^e^(*'0[x] 



for a sufficiently regular function x- 



Appendix F 

In the bosonized version (appendix A), we can write 
in (1761) 



H-d = J dsih 'ipl{s)dsil'i{s) - ipl{s)ds'4'r{s) 



(Fl) 
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with {■ijja{s),ij^l{s')} = 6abS{s — s') for a,b = l,r. To de- 
rive (1811) . we first introduce the auxiliary decomposition 

Idis) = iris) - ii{s), Qd{s) = Iri-S) + /^s) (F2) 

with the commutation rules 

[liis),ir{s')]^0, (F3) 

[ilis)Jiis')] = [Us')jr{s)]^ie^ds{s-S')/2TT. 

The last useful set of commutators is 

[Hd,Qd{s)\ = ihdsid{s), 
[HdJd{s)] ^ -ihdsQdis), 
[HdJr{s)] = ihdsiris), 
[HdJi{s)]^-ihdJiis). (F4) 



and 



Similarly to ([22|). equilibrium averages for the decoupled 
detector (A 0) are 

Tr/5d/d(a;) = 0, 
Tr/5d<9d(t^) = 2TmS{uj), 

TTpJi{a)ir{f3) = 0, (F5) 
TipJiia)ii{P) - TTpJr{l3)Ir{a) = 

= —6{a + l3){wd{a) + a), 

where n is average charge density and = 
w coth(fta;/fcBT'(;). We shall put n = assuming that 
the average charge is screened out. In the Heisenberg 
picture, we have 

{dt-ds)ii{s,t) = -dA{s)Q{t), 

The general solution is 

ds'X{s')Q{t + s~s'), 

ir{s,t) ^ ir{s-t) +ds ( ds'\{s')Q{t - S + s'). 

J —CO 

Hence, we have in the Heisenberg picture 

2/(s, t) = 2d, j ds'\{s')Q{t -\s^s'\)+ (F6) 

ir{s -t)- iri-S -t)~ii{t + S) +il{t~ S) 



ds X{s)Qdis,t) = / ds X{s)iiris~t)+ii{s + t)) 



-2 J ds I ds' \{s')X{s - s')Q{t - s). (F7) 



To find 5, we need the Keldysh generating functional,- 
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Tr/5 X 



(F8) 



^Jidt[xt){t)Q{t) + Sds{xt{s,t)i^{s,t)+xU^^t)h'^^M 

where — x/2 ± <t> a-nd S'[0, (j)] — 0. The functional is 
related to a quasiprobability in presence of external fields 



g[Ir,IuQ;<l)] = J Dxe- 



.s[x,0] X (F9) 

- f idt[xQ{t)Q{t)+f ds(Xr{sJ)Ir(s,t) + Xl{sJ)I,{s,t)}] 



One can show that Heisenberg equations (|F6[) and (|F7[) 
are satisfied for the quasiprobability and Ir, Ii and Q 
instead of the corresponding operators. It can be also 
shown that the backaction of the detector on the system 
can be simplified by the classical mapping Q — >■ Q, Qd — >■ 
Qd in Hamiltonian ((75)) used in (|F8)) . 

Hi{f) = ~2^n J ds\is){Qd{s,t)Q + Qd{s,t)Q{t))/e^. 

(FIO) 

The final generating functional reads 



.six] 



D<f>D£,DT] e 



jidt{mv{t)-x(i)m/'i) 



„-/rfu;[|x(i^)Pu.d(w)/327r^ + |0(w)|V2i«d(t^)] 

.S[~X*\-^,-2lT<j)*\+iTTfl] 



(Fll) 



with S defined by ^ and fi{t) ds\ * A(s)r;(i - s), 

it;(j(w) — cj coth(/iix)/A:sTc(). In the case of many tapes, 
([M)) follows from 



D(j)D£,Dri e 



/»<it(?(t),,(t)~X(t)Ee'?^<=(*)/2) 



X 



g-/<i^[|xMl'Ee"'cM/327r^ + Eel^cMlV2«'cM] x 
g^S[-(,-Y.a X*-^c,-27rX;c 0c*Ac+47ri)] (F12) 

POO 



with f]{t) = J2c /o°° ^^^c * K{s)ri{t - s) 
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